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Abstract 

A Brownian particle with diffusion coefficient D is confined to a bounded 
domain of volume V in R 3 by a reflecting boundary, except for a small 
absorbing window. The mean time to absorption diverges as the win- 
dow shrinks, thus rendering the calculation of the mean escape time a 
singular perturbation problem. We construct an asymptotic approxi- 
mation for the case of an elliptical window of large semi axis a <C V 1 ^ 3 

V 

and show that the mean escape time is Et ~ ^(e), where e is 

2irDa 

the eccentricity of the ellipse; and K(-) is the complete elliptic integral 

of the first kind. In the special case of a circular hole the result reduces 

V 

to Lord Rayleigh's formula Et ~ — — . which was derived by heuristic 

AaD 

considerations. For the special case of a spherical domain, we obtain 



V 

the asymptotic expansion Et 



a . R / a 
1 + — log — + O ( — 
R a \R 



This 



AaD 

problem is important in understanding the flow of ions in and out of 
narrow valves that control a wide range of biological and technological 
function. 



1 Introduction 

We consider the exit problem of a Brownian motion from a bounded domain, 
whose boundary is reflecting, except for a small absorbing window. The mean 
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first passage time to the absorbing window (MFPT), Et, is the solution of 
a mixed Neumann-Dirichlet boundary value problem (BVP) for the Poisson 
equation, known as the corner problem, which has singularity at the boundary 
of the hole EQ-jH]- The MFPT grows to infinity as the window size shrinks 
to zero, thus rendering its calculation a singular perturbation problem for the 
mixed BVP, which we call the narrow escape problem. The narrow escape 
problem has been considered in the literature in only a few special cases, be- 
ginning with Lord Rayleigh (in the context of acoustics), who found the flux 
through a small hole by using a result of Helmholtz 0. He stated [7] (p. 176) 
" Among different kinds of channels an important place must be assigned to 
those consisting of simple apertures in unlimited plane walls of infinitesimal 
thickness. In practical applications it is sufficient that a wall be very thin in 
proportion to the dimensions of the aperture, and approximately plane within 
a distance from aperture large in proportion to the same quantity. " More 
recently, Rayleigh's result was shown to fit the MFPT obtained from Brow- 
nian dynamics simulations [8 . Another result was presented in 0, where a 
two-dimensional narrow escape problem was considered and whose method is 
generalized here. A related problem is that of escape from a domain, whose 
boundary is reflecting, except for an absorbing sphere, disjoint from the reflect- 
ing part of the boundary ^TJ and references therein] . It differs from the narrow 
escape problem in that there is no singularity at the boundary and there is no 
boundary layer. The mixed boundary value problems of classical electrostatics 
(e.g., the electrified disk problem [12]), elasticity (punch problems), diffusion 
and conductance theory, hydrodynamics, and acoustics were solved, by and 
large, for special geometries by separation of variables. In axially symmetric 
geometries this method leads to a dual series or to integral equations that can 
be solved by special techniques The special case of asymptotic rep- 

resentation of the solution of the corner problem for small Dirichlet and large 
Neumann boundaries was not done for general domains. The first attempt in 
this direction seems to be j^j- The narrow escape problem does not seem to 
fall within the theory of large deviations jTH] • It is different from Kolmogorov's 
exit problem [F3\ of a diffusion process with small noise from an attractor of 
the drift (e.g., a stable equilibrium or limit cycle) in that the narrow escape 
problem has no large and small coefficients in the equation. The singularity of 
Kolmogorov's problem is the degeneration of a second order elliptic operator 
into a first order operator in the limit of small noise, whereas the singularity 
of the narrow escape problem is the degeneration of the mixed BVP to a Neu- 
mann BVP on the entire boundary. There exist precise asymptotic expansions 
of Et for Kolmogorov's exit problem, including error estimates (see, e.g., [3], 
jS]), which show that the MFPT grows exponentially with decreasing noise. In 
contrast, the narrow escape time grows algebraically rather than exponentially, 
as the window shrinks. Our first main result is a derivation of the leading or- 
der term in the expansion of the MFPT of a Brownian particle with diffusion 
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coefficient D, from a general domain of volume V to an elliptical hole of large 
semi axis a that is much smaller than V 1 ^ 3 , 



Et 



V 



K(e), 



(1.1) 



2nDa 



where e is the eccentricity of the ellipse, and K(-) is the complete elliptic 
integral of the first kind. In the special case of a circular hole (jl-ljl reduces to 



Equation (Jl.lj) shows that the MFPT depends on the shape of the hole, and not 
just on its area. This result was known to Lord Rayleigh |7], who considered 
the problem of the electrified disk (which he knew was equivalent to finding the 
flow of an incompressible fluid through a channel and to the problem of finding 
the conductance of the channel), who reduced the problem to that of solving 
an integral equation for the flux density through the hole. The solution of the 
integral equation, which goes back to Helmholtz [|J| and is discussed in [TH], is 
proportional to (a 2 — p 2 ) -1 / 2 in the circular case, where p is the distance from 
the center of the hole Note that equations and (jl.2|) are leading 

order approximations and do not contain an error estimate. We prove by 
using the singularity properties of Neumann's function for three-dimensional 
domains, in a manner similar to that used in jU] for two-dimensional problems. 
The leading order term is the solution of Helmholtz's integral equation [B]. 
Our second main result is a derivation of the second term and error estimate 
for a ball of radius R with a small circular hole of radius a in the boundary, 



Equation (jl.3j) contains both the second term in the asymptotic expansion of 
the MFPT and an error estimate. We use Collins' method EI] of solving 
dual series of equations and expand the resulting solutions for small e = a/R. 
The estimate of the error term, which turns out to be O (slogs), seems to be 
a new result. An error estimate for eq. (|l.l|) for a general domain is still an 
open problem. We conjecture that it is O (slogs), as is the case for the ball. 
If the absorbing window touches a singular point of the boundary, such as a 
corner or cusp, the singularity of the Neumann function changes and so do 
the asymptotic results. In three dimensions the class of isolated singularities 
of the boundary is much richer than in the plane, so the methods of |22] 
cannot be generalized in a straightforward manner to three dimensions. We 
postpone the investigation of the MFPT to windows at isolated singular points 
in three dimensions to a future paper. In Section |21 we derive a leading order 
approximation to the MFPT for a general domain with a general small window. 




(1.3) 
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The leading order term is expressed in terms of a solution to Helmholtz's 
integral equation, which is solved explicitly for an elliptical window. In Section 
El we obtain two terms in the asymptotic expansion of the MFPT from a ball 
with a circular window and an error estimate. Finally, we present a summary 
and list some applications in Section 0] This is the first paper in a series of 
three, the second of which considers the narrow escape problem from a bounded 
simply connected planar domain, and the third of which considers the narrow 
escape problem from a bounded domain with boundary with corners and cusps 
on a two-dimensional Riemannian manifold. 



2 General 3D bounded domain 

A Brownian particle diffuses freely in a bounded domain f2 C R 3 , whose bound- 
ary dQ is sufficiently smooth. The trajectory of the Brownian particle, denoted 
x(t), is reflected at the boundary, except for a small hole dQ a , where it is ab- 
sorbed. The reflecting part of the boundary is dQ r = dQ — dQ a . The lifetime 
of the particle in Q is the first passage time r of the Brownian particle from 
any point x G f2 to the absorbing boundary dQ a . The MFPT, 

v(x) = E[t\x(0) = x], 

is finite under quite general conditions [24J. As the size (e.g., the diameter) of 
the absorbing hole decreases to zero, but that of the domain remains finite, we 
assume that the MFPT increases indefinitely. A measure of smallness can be 
chosen as the ratio between the surface area of the absorbing boundary and 
that of the entire boundary, 



\dQ a \ 

\dn\ 



< i, 



(see, however, a pathological example in Appendix EJ) • The MFPT v(x) sat- 
isfies the mixed boundary value problem 



Av(x) = - — , for x G Q, (2.1) 



v{x) = 0, for x G dtt a , (2.2) 
dv(x) 



dn(x) 



0, for x G dVl r 



where D is the diffusion coefficient. According to our assumptions v(x) — ► oo 
as the size of the hole decreases to zero, e.g., as e — ► 0, except in a boundary 
layer near dVL a . Our purpose is to find an asymptotic approximation to v(x) 
in this limit. 
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2.1 The Neumann function and integral equations 

To calculate the MFPT v(x), we use the Neumann function N(x,£) (see jH], 
[11J), which is a solution of the boundary value problem 

&xN{x,£) = -6{x-£), for x,£etl, (2.3) 



dn{x) 



\dn\ 



for x e dn, £ e fi, 



and is defined up to an additive constant. The Neumann function has the form 



N{x,£) 



(2.4) 



An\x — £| 

where vs(x, £) is a regular harmonic function of x e and of £ G fi. Green's 
identity gives 



[JV(aj, ^)Au(aj) - v{x)AN{x,£)} dx 



on 



^( S ),«^-^(S))«« 

v(x(S))dS 



dn 



dn 



dQ\ 



on 



On the other hand, equations ()2.1|) and (|2.3|) imply that 



[JV(x,^)Av(x) -«(a;)AJV(x,£)] da; = 



N(x,g) dx, 



hence 



iV(a;,£)rfa3 



(2.5) 



v(x(S))dS. 



on 



Note that the second integral on the right hand side of eq. (|2.5jl is an additive 
constant. Setting 



C 



\dQ\ 



v{x{S))dS, 



(2.6) 



an 



we rewrite eq. (|2.5|) as 



-Yfar.O'/aM / N(xffl,t) ^*^ dS-C, (2.7) 



which is an integral representation of v(£). We define the boundary flux density 



S (S) = *!<=£>), ( ,s) 

choose £ G dQ a , and use the boundary condition (|2.2|) to obtain the equation 

= 1/ iV^r, £) da: + / iV(x(5), i)g(S) dS - C, (2.9) 
u Jn Jdn a 

for all £ G <9fi a . Equation (|2.9j) is an integral equation for g(S) and C. To 
construct an asymptotic approximation to the solution, we note that the first 
integral in equation ()2.9|) is a regular function of £ on the boundary. Indeed, 
due to symmetry of the Neumann function, we have from (|2.3j) 



At / N(x,£)dx = -1 for £eQ 
^ Jn 



[2.10) 



and 

d f „, ^ , m 



a»<eyo Ar( " ,t) ' to = -lsii for * em - (2 ' u) 

Equation (|2.1()jl and the boundary condition (|2.11j) are independent of the 
hole dQ a , so they define the integral as a regular function, up to an additive 
constant, also independent of dfl a . The assumption that for all x G Q, away 
from dfl a , the MFPT v(x) increases to infinity as the size of the hole decreases 
and eq. ([2.6)1 imply that C — > oo as as the size of the hole decreases to zero. 
This means that for £ G dQ a the second integral in eq. ()2.9)1 must also become 
infinite in this limit, because the first integral is independent of dQ a . Therefore, 
the integral equation (|2.9J) is to leading order 

/ N{x{S),t))g (S)dS = C for ^ G dVl a , (2.12) 

where go(S) is the first asymptotic approximation to g(S) and Co is the first 
approximation to the constant C. Furthermore, only the singular part of the 
Neumann function contributes to the leading order, so we obtain the integral 
equation 

1 f *<"> dS x = C , (2.13) 



2tt J d n a \x - y\ 

where Cq is a constant, which represents the first approximation to the mean 
first passage time (MFPT). Note that the singularity of the Neumann function 
at the boundary is twice as large as it is inside the domain, due to the con- 
tribution of the regular part (the "image charge" ) and therefore the factor — 

4:71 

of equation ()2.4|) was replaced by — . In general, the integral equation ()2.13|) 
has no explicit solution, and should be solved numerically. 
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2.2 Elliptic hole 

When the hole dfl a is an ellipse, the solution of the integral equation (j2.13|) is 
known jjj, [TH]. Specifically, assuming the ellipse is given by 

x^ iP' 

- + f^ = l, z=0, (b<a), 
a 2 ¥ 

the solution is 

9o(x) = 90 , (2.14) 

/ x y 

where go is a constant (to be determined below). The proof, originally given 
in 6j, is reproduced in Appendix El To determine the value of the constant 
go, we use the compatibility condition 

/ g {x)dS x = ^, (2.15) 
Jdn a u 

obtained from the integration of eq. ()2.1j) over Q. Using the value 



a 



g (x)dS x = I dx I pJ f ^ =2irabg (2.16) 

dn a J-a J-bsj\-K: ar y l 

V 

and the compatibility condition (|2.15j) . we obtain 

*-55*- (2 - 17) 

Hence, by equation ()B.5|) . the leading order approximation to C is 

C = -f dS x = %( e ), (2.18) 

2tt J dQa \x - y\ 2%Da 

where K(-) is the complete elliptic integral of the first kind, and e is the 
eccentricity of the ellipse, 

' IP 

(2-19) 

or 

In other words, the MFPT from a large cavity of volume \Q\ through a small 
elliptic hole is to leading order 

Er(a,b) ~ -^L K (e). (2.20) 
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71 

For example, in the case of a circular hole, we have e = and K(0) = — , so 
that 

Er{a,a)~M.=0(l), (2.21) 

provided 

I r> I 2/3 

0(1) for 



|<9fi| 

Equation ()2.21|) was used in [S], If the mouth of the channel is not circular, 
the MFPT is different. Equation (j2.21j) indicates that a Brownian particle that 
tries to leave the domain "sees" finer details in the geometry of the hole and the 
domain than just the quotient of the surface areas. The additional geometric 
features contained in the MFPT are illustrated by the two interesting limits 
e« 1, where the ellipse is almost circular, and 1 — e 1, where the ellipse 
is squeezed. In the case e <C 1, we use the expansion of the complete elliptic 
integral of the first kind 



(2 - 22) 

In the second limit 1 - e < 1, we find from the asymptotic behavior 



lim 

e-+l 



2 ° 1 - e 



(2.23) 



that 

I O I 1 f\ 

£r~^log , for 1 - e < 1. (2.24) 

4na 1 — e 

The area of the hole is given by 

S = nab = vraVl - e 2 , (2.25) 

or equivalently 

S l/2 

a = -7-, (2.26) 

and the MFPT has the asymptotic form 

& „«M! log JL, for (2.27, 
AVtTS 6 l-e' 
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3 Explicit computations for the sphere 



The analysis of Section [2] is not easily extended to the computation, or even 
merely the estimation of the next term in the asymptotic approximation of 
the MFPT. The explicit results for the particular case of escape from a ball 
through a small circular hole gives an idea of the order of magnitude of the 
second term and the error in the asymptotic expansion of the MFPT. If the 
domain Q is a ball, the method of [IB]- [IS], [20], and [21] can be used to 
obtain a full asymptotic expansion of the MFPT. We consider the motion of 
a Brownian particle inside a ball of radius R. The particle is reflected at the 
sphere, except for a small cap of radius a = eR and surface area 4nR 2 sin 2 -, 

where it exits the ball. We assume e < 1. The MFPT v(r,9,4>) satisfies the 
mixed boundary value problem for Poisson's equation in the ball [21], 

Av(r, 6,<f>) = -1, for r < R, < 9 < tt, < <p < 2tt, 



v(r,6,4>) 
dv(r, 9, 4>) 



Or 



0, for < 9 < e, < <p < 2vr, 
0, for e < 9 < tt, < <p < 2tt, 



(3.1] 



r=R 



The diffusion coefficient has been chosen to be D = 1. Due to the cylindrical 
symmetry of the problem, the solution is independent of the angle 0, that is, 
v(r, 9, <fi) = v(r, 9), so the system (j3.1j) can be written as 



Av(r 
v(r, 9) 

dv(r, 9) 



r=R 



dr 



-1, for r < R, < 9 < tt, 



0, for < 9 < e, 



0, for e < 9 < tt, 



r=R 



where the Laplacian is given by 
Av(r,8) 



1 d f odv* 
-^tt- r 2 — I - 



d 



The function f(r, 9) 



r 2 dr \ dr J r 2 sin 9 89 
R 2 -r 2 



Sin ^) 



6 

Af 



is the solution of the boundary value problem 
= —1, for r < R, 



f 



0. 



-R 



9 



In the decomposition v = u + f , the function u(r, 9) satisfies the mixed 
Dirichlet-Neumann boundary value problem for the Laplace equation 



Au(r, i 



u(r, 9) 
du(r, 9) 



0, for r < R, < 9 < tt, 



0, for < 9 < e, 



(3.2) 



Or 



r=R 



r=R 



R 



for e < 9 < 



7T. 



Separation of variables suggests that 

oo ^ 

u(r, 9) = ^a n P„(cos ( 



(3.3) 



ra=0 



where P n (cos^) are the Legendre polynomials, and the coefficients {a n } are to 
be determined from the boundary conditions 



u(r, 9) 
du(r, 9) 



dr 



r=R 



r=R 



COS( 



n=0 

oo 



y^na n P n \ 



0, < 9 < e 
R 2 



COS! 



e < 9 < 7r. 



(3.4) 



(3.5) 



n=l 



Equations ()3.4|) . (|3.5j) are dual series equations of the mixed boundary value 
problem at hand, and their solution results in the solution of the boundary 
value problem (|3.2j) . Dual series equations of the form 



^2a n P n {cos9) = 0, for < 9 < e, 



n=0 



^(2n + l)a n P n (cosfl) = G(6), for e < 9 < vr 



(3.6) 



(3.7) 



n=0 



are solved in O eqs.(5.5.12)-(5.5.14), (5.6.12)]. However, the dual series equa- 
tions ()3.6|) - (J3.7j) are different from equations ()3.4)l - ()3.5j) . The factor 2n+l that 
appears in equation ()3.7|) is replaced by n in equation ()3.5j) . What seems as a 
slight difference turns out to make our task much harder. The factor 2n + 1 
fits much more easily into the infinite sums (|3.6|) - (j3.7j) . because it is the nor- 
malization constant of the Legendre polynomials. 



3.1 Collins' method 

The solution of dual relations of the form ()3.5|) (see [121 (5.6.19)-(5.6.20)]) is 
discussed in [20], [21]. Specifically, assume that for given functions G{9) and 



10 



F{6) we have the representation 

oo 

^(l + F n )6 n T-™ n (cos^) = F(6), for < 9 < e, 



n=0 

oo 



J2( 2n + 2m + l)6 n T-™ n (cos^) = G(6), for e < 9 < vr, 



n=0 



where T m ™ n are Ferrer's associated Legendre polynomials |2Zj, [2H] and {P n } 
is a given series that is 0(n _1 ) as n ^ oo. Then for m = 0, we have 

oo 

^(l + #„)6 n P„(cos0) = for 0<£<£, (3.8) 

n=0 

oo 

5^(2n + l)6 n P n (cosfl) = G(9), for e < < vr. (3.9) 



n=0 



2n + 1 ___ 

Setting ao = &o, a n = b n , n > 1 in equations (I3.4I) - (I3.5|) results in 

2n 

oo 

5^(1 +H n )b n P n {cos9) = 0, for < 9 < e, (3.10) 

n=0 

°° 2P 2 

5^(2n + l)6 n P n (cos#) = + b 0} for £ < 9 < vr. (3.11) 

3 

n=0 

Equations (I3~TUJ) - (|3~TTJ) are equivalent to (j3~8 )l -([3~n j) with H Q = 0, H n = 
1 2P 2 

— , n > 1, P(#) = 0, and G{9) = h bo. Collins' method of solution 

2n 3 

consists in finding an integral equation for the function 

oo 

h(9) = J2(2n + l)b n P n (cos9), for < 9 < e, 

n=0 

so that 

if 6 i r 

b n — - / h(a)P n (cosa) sma da + - / G(a)P„(cos a) sin a da. 
2 Jo 2 J £ 

Substituting into equation ()3.8|) . with F{9) = 0, we find for < 9 < e that 

^ r-e 00 

= - / /j-(a) y^(l + P ra )P ra (cosa)P n (cos6 l ) sinacfa: 
2 -A) n=0 

+ 9 / G(a)^(l + P n )P n (cosa)P„,(cosfl)sinada. (3.12) 
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3.2 The asymptotic expansion 

To facilitate the calculations, we consider first the case H n = for all n. Then 
we will show that the leading order term obtained for this case is the same as 
that for the case H n ^ 0. In the latter case, we obtain the first correction to 
the leading order term and an estimate on the remaining error. 

3.2.1 The leading order term when H n = 

We will now sum the series (|3.12|) in the case H n = 0. First, we recall Mehler's 
integral representation for the Legendre polynomials |29j . 



where H(x) is the Heaviside unit step function. Then we obtain for u < 9 < 




(3.13) 



and the identity (THj 




(3.14) 



e < a, 




oo 



n=0 





(3.15) 



Similarly, 




oo 



n=0 




(3.16) 



Hence, 





(3.17) 
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Equation (J3.17)) means that the Abel transforms j3U] °f two functions are the 



same, so that 



h(a) sin a da 
u \J cos u — cos a 



w G(a) sin a da 
\J cos u — cos a ' 



(3.18) 



because the Abel transform is uniquely invertible. Equastion ()3.18|) is an Abel- 
type integral equation, whose solution is given by 



h{ff)f*nO = ~ 



sinudu f n G(a) sin a da 



vr d9 J d y/cos 9 — cos u J e V cos u — cos a 



or 

where 
and 



h(0) 



2 d 



H(u) sinudu 



G(u, e 



sin 9 d9 J e y/ cos 9 — cos u 
H{u) = -G(u,e), 

1 r G(9)sin9d9 



2vr J £ \J cos u — cos 9 



(3.19) 

(3.20) 
(3.21) 
(3.22) 



2R 2 



The dual integral equations ()3.1U|) - ([3.1ip define G{9) = h bo, so that 

3 



1 

2^ 



2R 2 



+ b 



sin 9 d9 



2R 2 
2R 2 



\J cos ip — cos 9 

+ b(A — a/ cos ■?/> — cos 9 
J * 

+ 6 )I(V2c<4- 



(3.23) 



cos -0 — cos , for ^ < 



In particular, setting n = in equation ()3.12|) and using equation ()3.20|) . gives 



b 



h(a) sin a da + 



1 r (2R 2 



b I sin a da (3.24) 



V2 [ H(tjj) cos - dip + 



2i? 2 



+ 6o cos 
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Integrating equation (|3.23|) . we obtain 

V2 [ G(tp,e) cos^ dtp = (3.25) 
Jo 2 

[2R 2 L \ V2 T / r- i> , ■ \ V j , 

= h bo — / v 2 cos y cos-^ — cose cos — dip 

V 3 / n Jo V 2 / 2 

2i? 2 , £ 
T + S ^ ^ t \ 4 /- Sm 2 S 2 rf S 



7T \ 3 / 7T 







sin s 

2 



2R 2 

(e + sine) — — — h o sin -. 



n ■ -— , . 3 ■ o , — 2 - 



Combining equations ()3.24j) and (j3.25j) gives 

2i? 2 / 7T \ 2R 2 / n _\ \n 



a 



3 \e + sine / 3 \2e J 4a V Vi? 

where |Q| = — - — is the volume of the ball, and a = Re is the radius of the 
hole. 

3.2.2 The case H n ^ 

The asymptotic expression (|3.26|) for bo, was derived under the simplifying 
assumption that H n = 0. However, we are interested in the value of bo which 
is produced by the solution of the dual series equations ()3.1U|) - (|3.11|I . where 

H n = — . We sum the series (I3.12|) by the identities 

(cosa)P n (cos6) sin a da 

2 J o n= o 

1 r,/ v 72 r c ° s ( n + k) v dv V 2 f 6 cos(n+ \)udu 

= - h(a) y H n — / —^^^=^^^ — / sin a da 

2 Jo ~£ 71 Jo vcos^-cos" n Jo Vcosw-cosfl 



1 r r 

— / h(a) sin a da I 
2 tt Jo Jo 



dv f e K(u, v) du 



V cos v — cos a Jo v cos w — cos 6 



du f e TW . , f e h(a) sin a da 

' K(u,v)dv / ; v ' =, (3.27) 



2 7T Jo V COS M — COS 9 Jo Jv V cos v ~ cos a 
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where 



K(u, v) 



— H n cos ( n -\ — )ucos(n-\ — )v (3.28) 



n=0 

y log2 



+ 



2n 



cos ^(v — u) 



2tt 

V + U — 7T 1 



log 2 



sin -(f + w, 



sin — (u — u. 
2 V 



47T 



— — ^- I 

sin - Ou + u) H sin -(t> — it) 

2 47r 2 



Similarly, 



^^if n P n (cos a:)P n (cos 0) sinarfa 



ra=0 



1 /" r f c 
— / G[a) sin a da I 
2tt J £ Jo 



du 



27T Jo V cos u — cos 6> 



v) du 

y/ cos t> — cos a J \/ cos u — COS # 



sin a da 



v/cosf — cos a; 



. (3.29) 



Substituting equations JHHSl, IpTTBjl . ()T77jl . and (jHIZHD into equation (j3~T^|l 
yields 



1 

= — 



du 



27T Jo V COS M — cos 9 Ju \/cosu — cos a 



h(a) sin a da 



+ - 



1 



du 



27T Jo V cos u — COS 6> Jo 



K(u, v) dv 



h(a) sin a da 
y/cos v — cos a 



1 

+ 7^ 



77 G(a)sinao?a 

27T 7o V COS U — COS J E \/cOSU — COS 



+ - 



27T Jo V cos u — COS 6* 



G(a) sin a da 



K(u, v) dv 
y/cos v — cos a ' 



which is again an Abel-type integral equation. Inverting the Abel transform 
|3*Uj . we obtain 







h(a) sin a da 1 
27r J u y/cosu — cos a 27r 



1 Z" 77 G(a)sina<ia 1 
27r J e ^/ cos u — cos a 27r 



K(u, v) dv 



h(a) sin a da 



\/ cos v — cos a 



(3.30) 



G(a) sin a da 



a K(u,v)dv 
V cos t> — cos a 
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Setting 



H{u) = ±- [ £ h M^t , (3.31) 
we invert the Abel transform ()3.31|) to obtain 

h[9) = -~ r *™ UH{u)dU ■ (3.32) 
sin & c/6 1 y e Vcos — cos u 

Writing 

J(«) = f/"(«)+G(«,e), (3.33) 

equation (|3.3U|) becomes 



J(«)+ / iif(tt,t;)J(7;)dv = M(«), (3.34) 
Jo 

where the free term M(u) is given by 

K(u,v)G(v,v)dv. (3.35) 
Equation ()3.34|) is a Fredholm integral equation for J. 

3.2.3 The second term and the remaining error: L 2 estimates 
Equations flUSH , flH251) , and (j3~33|) give that 

2i? 2 2R 2 tt v^tt f u 

6 + ^ = — — — + — / J{u)cos-du : (3.36) 

3 S e + sme e + sin e J 2 

where J is the solution of the Fredholm equation ()3.34|) . In this section we 
show that 

J(tt) cos — au = I o H — — I £ log — h C/(e / 



e + sin eJo 2 \ 3 / \ e 

therefore the last term in eq. (j3.36|) should be considered a small correction 

R 2 tt 

to the leading order term , obtained in Section Hi 2. 21 This confirms the 

3 e 

intuitive results of jH], |2Sj and gives an estimate on the error term. Due to 
the logarithmic singularity of the function K(u, v) (see (13.28)) ) the operator K, 
defined by 

Kf{u)= f K(u,v)f(v)dv, (3.37) 
Jo 

maps L 2 [0,e] into L 2 [0,e]. In Appendix 1X1 we derive the estimate 

||K|| 2 <^£log-, (3.38) 

for e < 1. Better estimates can be found; however we settle for this rough 
estimate that suffices for our present purpose. 
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3.2.4 Estimate of \\j \\ 2 

In terms of the operator K, equation ()3.34j) can be written as 

J = M- KJ. (3.39) 

The triangle inequality yields 

II Jh < \\ M h + \\KJh < \\M\\ 2 + \\K\\ 2 \\ J\\ 2 , (3.40) 
which together with the estimate (|3.38|) gives 

< : ^ M }tu < ( 1+slog-l ||M|| 2 for £<1. (3.41) 



l-\\Kh ~ V °£ 



3.2.5 Estimate of ||M|| 2 

We proceed to find an estimation for ||M||2. First, we prove that the kernel 
satisfies the identity 

r v 

l K(u, v) cos - dv = 0, for all u. (3.42) 
Jo 2 

Indeed, by changing the order of summation and integration, we obtain 

i^i / i\ r ( i\ « 

— > — cos n H — I u I cos \ n H — v cos - dv 



v 

K(u, v) cos - dv 



— — COS ( n + — ] ii I fcosf// -4 1 )r -f cos //r ) c/r 



2vr ^ n V 2 

ra=l 



= 0. (3.43) 
Equations (|3~2l?J) . (JSIS3), and (l3~4H imply that 

M(m) = — + &ol [ £ K(u,v) cos- dv. (3.44) 

V 3 /Jo 2 

The estimate ()3.38|) gives 

y/2 (2R 2 ,\ u „,, r- VTE f2R 2 , \ ,,„, 1 
||M|| 2 < ^- f — + b j \\K\\ 2y /i < ^- f — + 6 J ^ 3/2 log -• (3.45) 

Combining the estimates (j3.41|) and ()3.45|) . we obtain for e <C 1 

4 /2i? 2 \ 1 /2i? 2 , \ _ , /2l 

||J|| 2 <^(-g- + &o) £ 3/2 log-= ( — + M 0( £ 3 / 2 log £ ). 



(3.46) 
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3.2.6 The second term and error estimate 

The Cauchy-Schwartz inequality implies that 



71 



u 



J{u) cos — du 



< 



e + sin e 

for e C 1, which together with ()3.36|) gives 
bo = ^(l + 0(6loge)) 



2R2 u \ i 1 

— - + 6q e log 



(3.47) 



3c 



|0| 
4a 



;i + 0(eloge)). 



(3.48) 



To obtain the explicit expression for the term 0(e logs), we write the Fredholm 
integral equation (|3.34j) as 

(I + K)J = M. (3.49) 

The estimate ()3.38|) implies that \\K\\2 < 1 for sufficiently small e, hence 

J = M + O (\\K\\ 2 \\M\\ 2 ) . (3.50) 

Thus, using equation (|3.44|) and the estimates ()3.38|) and (|3.45|) . we write the 
last term in equation (|3.36|) as 



u 



u 



J{u) cos- du= / M(u) cos - du + O (e\\K\\ 2 \\M\\2 



(3.51) 



V2 ( 2R" 



7T V 



I / K(u, v) cos — cos — du dv + O {s^ log 2 
Jo Jo 2 2 



Equation ()3.28j) gives the double integral as 



o 



u v . . 1 2 1 . 2 n 



K(u, v) cos — cos -dudv = —e log — \- 0(e 
2 2 71 e 



hence 



y/2- 



71 



e + sin e 



u 



J(u) cos — du = I bo + 



2tf 



e log - + 0(e) 



Now it follows from equation ()3.36|) that 



bo 



M 

4o" 



1 + £-log-J- +0(e) 



(3.52) 
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3.3 The MFPT 



Using the explicit expression (|3.52|) . we obtain the MFPT from the center of 
the ball as 



u 



r=0 



r=0 



R 2 , R 2 

+ T = bo + Y = ^ 



l + elog- + 0(£) 

e 



(3.53) 



This is also the averaged MFPT for a uniform initial distribution, 



Et 



2- 



R 



sin6d8 / v(r,8)r dr 



|0| 
4a 



1 + slog - + 0(e) 



4 Summary and applications 

The narrow escape problem for a Brownian particle leads to a singular pertur- 
bation problem for a mixed Dirichlet-Neumann (corner) problem with large 
Neumann part and small Dirichlet part of the boundary. The corner problem, 
that arises in classical electrostatics (e.g., the electrified disk), elasticity (punch 
problems), diffusion and conductance theory, hydrodynamics, acoustics, and 
more recently in molecular biophysics, was solved hitherto mainly for special 
geometries. In this paper, we have constructed a leading order asymptotic ap- 
proximation to the MFPT in the narrow escape problem for a general smooth 
domain and have derived a second term and an error estimate for the case of 
a sphere. Our derivation makes Lord Rayleigh's qualitative observation into 
a quantitative one. Our leading order analysis of the general case uses the 
singularity property of the Neumann function for a general domain in R 3 . The 
special case of the sphere is analyzed by a method developed by Collins and 
yields a better result. A different approach to the calculation of the MFPT 
would be to use singular perturbation techniques. The vanishing escape time 
at the boundary would then be matched to the large outer escape time of order 
e^ 1 by constructing a boundary layer near the boundary. The analysis of the 
MFPT to a small window at an isolated singular point of the boundary is post- 
poned to a future paper. Brownian motion through narrow regions controls 
flow in many non-equilibrium systems, from fluidic valves to transistors and 
ion channels, the protein valves of biological membranes [31] . Indeed, one can 
view an ion channel as the ultimate nanovalve — nearly picovalve — in which 
macroscopic flows are controlled with atomic resolution. In this context, the 
narrow escape problem appeared in the calculation of the equilibration time of 
diffusion between two chambers connected by a capillary [23] • The equilibra- 
tion time is the reciprocal of the first eigenvalue of the Neumann problem in 
this domain, which depends on the MFPT of a Brownian motion in each cham- 
ber to the narrow connecting channel. The first eigenfunction is constructed 
by piecing together the eigenfunctions of the narrow escape problem in each 
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chamber and in the channel so that the function and the flux are continuous 
across the connecting interfaces. It was assumed in [23] that the flux profile 
in the connecting hole was uniform. The structure of the flux profile, which 
is proportional to (a 2 — p 2 ) -1 / 2 , has been observed by Rayleigh in 1877 [JJ. 
Rayleigh first assumed a radially uniform profile of flux and then refined the 
profile of flux going through the channel, allowing it to vary with the radial 
distance from the center of the cross section of the channel, so as to minimize 
the kinetic energy. A calculation of the equilibration time was carried out in 
[32*] by solving the same problem, and gave a result that differs from that of [JJ , 
which was obtained by heuristic means, by less than two percent. A different 
approximation, based on the Fourier-Bessel representation in the pore, was de- 
rived in [T<5] • Another application of the narrow escape problem concerns ionic 
channels |31j . and particularly particle simulations of the permeation process 
[33]- [37] that capture much more detail than continuum models. Up to now, 
computer simulations are inefficient because an ion takes so long even to enter 
a channel and then so many of the ions return from where they came. From the 
present analysis, it becomes clear why ions take so long to enter the channel. 
According to (jl.2|) the mean time between arrival of ions at the channel is 

Et 1 . . 

where N is the number of ions in the simulation and C is their concentration. 
A coarse estimate of f at the biological concentration of O.lMolar, channel 
radius a = 20A, diffusion coefficient D = 1.5 X lCT 9 m 2 /sec is r « lnsec. In 
a Brownian dynamics simulation of ions in solution with time step which is 
10 times the relaxation time of the Langevin equation to the Smoluchowski 
(diffusion) equation at least 1000 simulation steps are needed on the average 
for the first ion to arrive at the channel. It should be taken into account that 
most of the ions that arrive at the channel do not cross it [SB] . 

The narrow escape problem comes up in problems of the escape from a 
domain composed of a big subdomain with a small hole, connected to a thin 
cylinder (or cylinders) of length L. If ions that enter the cylinder do not return 
to the big subdomain, the MFPT to the far end of the cylinder is the sum of 
the MFPT to the small hole and the MFPT to the far end of the narrow 
cylinder. The latter can be approximated by a one-dimensional problem with 
one reflecting and one absorbing endpoint. If the domain has a volume V, the 
approximate expression for the MFPT is 

V L 2 

Et « -— + — . (4.2) 
AeD 2D v ; 

This method can be extended to a domain composed of many big subdomains 
with small holes connected by narrow cylinders. The case of one sphere of 
4irR 3 

volume V = , with a small opening of size e connected to a thin cylinder 
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of length L is relevant in biological micro-structures, such as dendritic spines 
in neurobiology. Indeed, the mean time for calcium ion to diffuse from the 
spine head to the parent dendrite through the neck controls the spine-dendrite 
coupling [35]. This coupling is involved in the induction of processes such 
as synaptic plasticity [10]. Formula ()4.2|) is useful for the interpretation of 
experiments and for the confirmation of the diffusive motion of ions from the 
spine head to the dendrite. 

Another significant application of the narrow escape formula is to provide 
a new definition of the forward binding rate constant in micro-domains [TU] . 
Indeed, the forward chemical constant is really the flux of particles to a given 
portion of the boundary, depending on the substrate location. Up to now, the 
forward binding rate was computed using the Smoluchowski formula, which 
corresponds to the absorption flux of particles in a given sphere immersed in 
an infinite medium. The formula applies when many particles are involved. 
But to model chemical reactions in micro-structures, where a bounded domain 
contains only a few particles that bind to a given number of binding sites, the 
forward binding rate, 

^forward ~~ ~I> 
has to be computed with f given in eq. ()4.1|) . 



A Estimate of II XI 



A.l Estimate of the kernel 



A rough estimate of the kernel, for < u, v < e, is obtained from equation 
(EPSfr as 

1 \ 2 
sin — (v + u) 



K 2 (u,v) < — cos-(f + u) ( log 2 



5 1 
+ —, ^cos-(v - u) ( log 2 

47T 2 2 



sin - ( v — u) 



Furthermore, 



£ 1 
cos — (v + u) I log 2 



sin - (v + u) 



du 



2 sin \{v+e) 



(logo;) 2 dx 



< 2 j sin -(v + e) — sin -v j j log 2 



1 

sin —v 
2 



< e cos —v ( log 2 



1 

sin —v 
2 



and 



r 1 

/ £ COS - 

Jo 2 



v ( log 2 sin —v I dv = e 



2 sin ie 



(\ogx) z dx < 2e 2 log 2 £:. 
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Similarly, 



£ 1 

COS ~(v — U) ( loj 

n„i„ 1.. 



2 sin — (v — u) 
2 V 



eft* 



2 r2sm±(e-u) 

(log x) dx + I (log x) efx < 

o Jo 

2u log 2 w + 2(e — u) log 2 (£ — «)• 



It follows that 



/ (2u log 2 u + 2(e - u) log 2 (e - «)) < As 2 log 2 e, 
Jo 

because ulogw is an increasing function in the interval < u < e~ 2 . Alto- 
gether, we obtain 



\K\\o < 



30 . 1 . _ 2 
e log - for e < e , 

27T £ 



(A.1) 



which is (jUSHD - 



B Elliptic hole 

We present here, for completeness, Lure's [TH] solution to the integral equation 
(I2.13J1 in the elliptic hole case. We define for y = (x, y) 

Hy) = l-- 2 ~i (b<a) 
and introduce polar coordinates in the ellipse dQ a 

x — y + (p cos 9, p srn 9) ) 

with origin at the point y. The integral in eq. (j2.1Hj) takes the form 

* M9) ~go dp 



9o{x) 

an a \ x — y\ 



dS T 



2n 

de 

JO 



(B.l) 



where po(9) denotes the distance between y and the boundary of the ellipse in 
the direction 9. Expanding L(x) in powers of p, we find that 



L(x) = 1 - 



(x + pcos6*) 2 (y + psim 
a 2 b 2 



L{y) - 24> lP - 2 p 2 , (B.2) 
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x cos 9 y sin cos 2 6 sin 2 6 
where <pi = 1 — — and q>i = — 1 — — . Solving the quadratic 



b 2 



b 2 



equation (|B.2|) for p, taking the positive root, we obtain 

p(x) = i + [0 2 + 2 (L(y) - L(x))] 1/2 



therefore, for fixed y and 9, 
dp(x) = 



1 



dL{x) 



2[0 2 + 2 (L(y) - L(x))) 1/2 ' 
and the integral takes the form 



(B-3) 



(B.4) 



dS r 



2- 



de 



dL(x) 



9o 



2- 



dB 



L(y) i 

^W+<t> 2 {L{y)-L{x)f /2 ^L{ 
g dz 



Substituting s 



,(x s 



and setting ip 



dS T 



an a \x-y 

2tt 



2tt 



de 



go 



- — -, we find that 
ds 



-i 



2 IT 



d6^2, rctan J t±JL 
1 — s 



9o 



2\/02 Jo \P$ + Sy/1 - S 
1 







^7r — 2 arctan \ dd 



( 



277 



<?o dO 



cos 2 6 1 sin 2 6 1 



x cos 6* y sin # \ 



+ 



7r — 2 arctan 



6 2 



cos 2 sin 2 6 



+ 



0^ \ V 6 2 

The arctan term changes sign when 9 is replaced by 8 + tt, therefore its integral 
vanishes, and we remain with 

»27T 



2 



27r6p 



d0 



cos 2 8 sin 2 6 1 



b 2 



de 



i) 2 



■ sin 2 6 



2nbg K(e), 



(B.5) 
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where K(-) is the complete elliptic integral of the first kind, and e is the 
eccentricity of the ellipse 



b 2 

- -, (a > b). (B.6) 

a z 

We note that the integral (jB.5|) is independent of y, so we conclude that (j2.14j) 
is the solution of the integral equation (|2.13j) . 



C A pathological example 

We have derived an integral equation for the leading order terms of the flux and 
the MFPT in the case where the MFPT increases indefinitely as the relative 
area of the hole decreases to zero. However, the MFPT does not necessarily 
increase to infinity as the relative area of the hole decreases to zero. This 
is illustrated by the following example. Consider a cylinder of length L and 
radius a. The boundary of the cylinder is reflecting, except for one of its 
bases (at z = 0, say), which is absorbing. The MFPT problem becomes one 
dimensional and its solution is 

v(z) = Lz-£. (C.l) 

Here there is neither a boundary layer nor a constant outer solution; the MFPT 
grows gradually with z. The MFPT, averaged against a uniform initial dis- 

L 2 

tribution in the cylinder, is Er = — and is independent of a, that is, the 
assumption that the MFPT becomes infinite is violated. 
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